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Abstract. The energy spectrum of nuclear recoils in Weakly Interacting Massive 
Particle (WIMP) direct detection experiments depends on the underlying WIMP mass. 
We study how the accuracy with which the WIMP mass could be determined by a single 
direct detection experiment depends on the detector configuration and the WIMP 
properties. We investigate the effects of varying the underlying WIMP mass and cross- 
section, the detector target nucleus, exposure, energy threshold and maximum energy, 
the local circular speed and the background event rate and spectrum. The number 
of events observed is directly proportional to both the exposure and the cross-section, 
therefore these quantities have the greatest bearing on the accuracy of the WIMP 
mass determination. The relative capabilities of different detectors to determine the 
WIMP mass depend not only on the WIMP and target masses, but also on their energy 
thresholds. The WIMP and target mass dependence of the characteristic energy scale 
of the recoil spectrum suggests that heavy targets will be able to measure the mass 
of a heavy WIMP more accurately. We find, however, that the rapid decrease of the 
nuclear form factor with increasing momentum transfer which occurs for heavy nuclei 
means that this is in fact not the case. Uncertainty in the local circular speed and 
non-negligible background would both lead to systematic errors in the WIMP mass 
determination. For deviations of ±20kms~^ in the underlying value of the circular 
speed the systematic error is of order 10%, increasing with increasing WIMP mass. 
This error can be reduced by also fitting for the circular speed. With a single detector it 
will be difficult to disentangle a WIMP signal (and the WIMP mass) from background if 
the background spectrum has a similar shape to the WIMP spectrum (i.e. exponential 
background, or flat background and a heavy WIMP). 
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1. Introduction 



Cosmological observations indicate that the majority of the matter in the Universe is 
dark and non-baryonic (e.g. Ref. [1]). Weakly Interacting Massive Particles (WIMPs) 
are one of the leading cold dark matter candidates, and supersymmetry provides a 
concrete, well-motivated WIMP candidate in the form of the lightest neutralino (e.g. 
Ref. [21 [3] ) . The direct detection of WIMPs in the lab [1] (see Ref. [5] for a review of 
current and future experiments) would not only directly confirm the existence of dark 
matter but would also allow us to probe the WIMP properties, as the shape of the 
differential event rate depends on the WIMP mass P [TJ El El [IDl [H]. Constraints 
on, or measurements of, the WIMP mass and elastic scattering cross-section will 
be complementary to the information derived from collider and indirect detection 
experiments |12j . 

In paper I [10], see also Refs. [9l HI], we examined the accuracy with which a 
future SuperCDMS [I3] like direct detection experiment would be able to measure 
the WIMP mass, given a positive detection. In this paper we revisit that analysis, 
studying in more detail the dependence of WIMP mass limits on the detector capabilities 
(including threshold energy, exposure, maximum energy and target nucleus), the WIMP 
properties (mass and cross-section), the local circular speed and the effects of non-zero 
backgrounds. We outline the calculation of the differential event rate and the Monte 
Carlo simulations in Sec. [2] (see paper I [10] for further details), present the results in 
Sec. [3] and conclude with discussion in Sec. HI 



2. Method 



2.1. Event rate calculation 

The differential event rate, assuming spin-independent coupling, is given by (see 
e.g. [21 E]): 

^(E) = -^A^F^E^E) , (1) 

where p-^ is the local WIMP density, the WIMP scattering cross section on the 
proton, /ip^ = (mpm^)/(mp + m^) the WIMP-proton reduced mass, A and F{E) the 
mass number and form factor of the target nucleus respectively and E is the recoil 
energy of the detector nucleus. We use the Helm form factor [14j. The dependence on 
the WIMP velocity distribution is encoded in J-'{E), which is defined as 

HE) = { r ^^dv) , (2) 



where f^{v,t) is the (time dependent) WIMP speed distribution in the rest frame of 
the detector, normalized to unity and (..) denotes time averaging. The WIMP speed 
distribution is calculated from the velocity distribution in the rest frame of the Galaxy, 
/*^(v), via Galilean transformation: v ^ v = v + v^(t) where v^(t) is the Earth's 
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velocity with respect to the Galactic rest frame [B]. The lower limit of the integral, fmin, 
is the minimum WIMP speed that can cause a recoil of energy E: 

1/2 

(3) 



EmA 



where is the atomic mass of the detector nuclei and /xax the WIMP-nucleon reduced 
mass. We use the 'standard halo model', an isotropic isothermal sphere, for which the 
local WIMP velocity distribution, in the Galactic rest frame, is Maxwellian 

/G(v) = N [exp (-|v| Vt;2) - exp {-viJvl)] |v| < v,,, , (4) 

/''(v) = |v|>t;esc, (5) 

where is a normalization factor and Vc = 220±20 km s"^ [15j and Vgsc ~ 540 km s~^ [16] 
are the local circular and escape speeds respectively. If the ultra-local WIMP 
distribution is smooth, then the uncertainties in the detailed shape of the local velocity 
distribution lead to relatively small changes in the shape of the differential event 
rate [T7l|18]. Consequently there is a relatively small, [C(5%)], systematic uncertainty 
in the WIMP mass [TO]. We caution that the assumption of a smooth ultra-local WIMP 
distribution may, however, not be valid on the sub milli-pc scales probed by direct 
detection experiments (e.g. Ref. [19] but see also Ref. [20] for arguments that the ultra- 
local WIMP distribution consists of a large number of streams, and is hence effectively 
smooth). 

As shown by Lewin and Smith p], see also Paper I [10], the differential event rate 
can, to a reasonable approximation, be written as 

The event rate in the E — >• OkeV limit, (di?/di?)o, and i?R, the characteristic energy 
scale, are given by 

§) =coy^^^ (7) 

and 

E. = c,JJ^, (8) 

respectively, where Cq and ce^^ are constants of order unity which are required when the 
Earth's velocity and the Galactic escape speed are taken into account and are determined 
by fitting to the energy spectrum calculated using the full expression, eq. ([1]). The exact 
values of these constants depend on the target nucleus, the energy threshold and the 
Galactic escape speed. For a Ge detector with energy threshold Eth = OkeV, cq ~ 0.78 
and ce^ ~ 1.72, with a weak dependence on the WIMP mass [lOj. For the majority of 
our calculations we will use the accurate expression, eq. ([I]), however in Sec. 13. 31 where we 
consider varying Vc we will use the fitting function, eq. ([6]), as it is not computationally 
feasible to carry out the full calculation in the likelihood analysis in this case. 
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Figure 1. From top panel to bottom: the expected number of events, A'oxp, for a 
£ = 3 X 10'^ kg day exposure (with cTp = 10~^pb), the characteristic energy scale, -Er, 
and the variation of the characteristic energy scale with mass, dEfi/dm^, as a function 
of WIMP mass for a Ge detector (solid line) and a Xe detector (dashed) . For A'exp the 
the energy thresholds of the current CDMS II [21] and XenonlO [22] experiments have 
been used: i?th = 10 keV and 4.5 keV for Ge and Xe respectively. 



The characteristic energy, i^R, depends on the WIMP mass, m^^, and the mass of 
the target nuclei, mp^. For WIMPs which are hght compared with the target nuclei, 

<C mA, -Er oc m^/mp^, while for heavy WIMPs, ^ uia, -Er ~ const. In words, 
for light WIMPs the energy spectrum is strongly dependent on the WIMP mass while 
for heavy WIMPs the dependence on the WIMP mass is far weaker. Consequently 
it will be easier to measure the mass of light (compared with the target nuclei) than 
heavy WIMPs. Since the experiments which currently have the greatest sensitivity are 
composed of Ge and Xe (CDMS II [21] and XenonlO [22] respectively) we focus on these 
targets. Fig. [1] shows the dependence of the characteristic energy, E^, and dE^/dm^ on 
the WIMP mass for Ge and Xe. For < (>) ~ 50 GeV varies more strongly with 

for Ge (Xe), reflecting the asymptotic WIMP mass dependences of the expression for 
i?R. This suggests that light (heavy) target nuclei will be better suited to determining 
the mass of light (heavy) WIMPs (see however Sec. 13.21) . The detector energy threshold 
will also come into play, in particular for small exposures, as the expected number of 
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events depends on the energy threshold. The expected number of events, A^'exp, for an 
exposure of £^ = 3 x 10^ kg day is also shown in fig. [1] as a function of WIMP mass for Ge 
and Xe detectors with Eth = lOkeV [211 113] and 4.5 keV [22] respectively. With these 
thresholds A^exp is larger (by a factor of order ~ 2 for ~ 50 GeV) for Xe than for 
Ge. This indicates that the relative capabilities of detectors to determine the WIMP 
mass will depend not only on the WIMP and target masses, but also on their energy 
thresholds. 



2.2. Monte Carlo simulations 

We use Monte Carlo simulations to examine, for a range of detector configurations and 
input WIMP masses, how well the WIMP mass could be determined from the energies 
of observed WIMP nuclear recoil events. 

We estimate the WIMP mass and cross-section by maximizing the extended 
likelihood function (which takes into account the fact that the number of events observed 
in a given experiment is not fixed), e.g. Ref. |23j : 



XN^^pt exp (-A). 



L = - v ^'^ utrfiEd . (9) 

Here Ne^^t is the number of events observed, Ei{i = 1, A^expt) are the energies 
of the events observed, f{E) is the normalized differential event rate and A = 
£ J^^(di?/di?) dE is the mean number of events where £ is the detector exposure (which 
has dimensions of mass times time) and ii^th is the threshold energy. We calculate the 
probability distribution of the maximum likelihood estimators of the WIMP mass and 
cross-section, for each detector configuration and input WIMP mass, by simulating 10^ 
experiments. We first calculate the expected number of events, Ajn, from the input 
energy spectrum. The actual number of events for a given experiment, A^^expt, is drawn 
from a Poisson distribution with mean Ain. We Monte Carlo generate Nexpt events 
from the input energy spectrum, from which the maximum likelihood mass and cross- 
section for that experiment are calculated. Finally we find the (two-sided) 68% and 95% 
confidence limits on the WIMP mass from the maximum likelihood masses. 



3. Results 



In Sec. 13.11 we investigate the mass limits for a a SuperCDMS like Ge detector [13] and 
their dependence on the detector energy threshold, maximum energy and exposure, and 
the WIMP cross- sect ion. In Sec. 13.21 we compare the Ge mass limits with those for a 
Xe detector, before examining the effects of uncertainties in the local circular speed and 
non-negligible background in Sees. 13.31 and 13.41 respectively. 



3.1. Germanium 

We begin, as in paper I \IQ\, by looking at a SuperCDMS like detector [13j, composed 
of Ge with a nuclear recoil energy threshold of i?th = lOkeV. We assume that the 
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Figure 2. Limits on the WIMP mass, m^™, as a function of the input WIMP 
mass, m™, for the benchmark Ge detector (i?th = lOkeV, perfect energy resolution, 
no upper hmit on energy of events detected and zero background) for exposures 
£ = 3 X 10"^, 3 X 10'' and 3 x 10^ kg day and input cross-section dp = 10"^ pb. The dot- 
dashed hnc is the input mass and the sohd (dotted) hncs are the 95% (68%) confidence 
limits. 



background event rate is negligible, as is expected for this experiment located at 
SNOLab [13], and that the energy resolution is perfect 0. For simplicity we assume that 
the nuclear recoil detection efficiency is independent of energy. The energy dependence 
of the efficiency of the current CDMS II experiment is relatively small (it increases 
from ~ 0.22 at E = Eth = lOkeV to ~ 0.30 at E = 15keV and then remains roughly 
constant). For further discussion of these assumptions see Ref. [TO] . 

We consider fiducial efficiency weighted exposures [§| £^ = 3 x 10^, 3 x 10^ and 
3 X 10^ kg day which correspond, roughly, to a detector with mass equal to that of 
the 3 proposed phases of SuperCDMS taking data for a year with a ~ 50% detection 
efficiency ]]]. 

I Gaussian energy resolution, with full width at half maximum of order 1 keV [13], does not affect the 
WIMP parameters extracted from the energy spectrum [lOj . 

§ For brevity we subsequently refer to this as simply the exposure. 

II The 50% detection efficiency was chosen based on Ref [24]. The more recent CDMS II analysis [21] 
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Figure 3. As fig. [2 but with the fractional deviation of the WIMP mass hmits from 
the input mass, (m'^™ — m™)/ra™, plotted. 

We use fiducial values for the detector energy threshold and WIMP-proton cross- 
section of E^\^ = lOkeV and (Xp = 10^^ pb but later, in Sees. 13.1.11 and 13.1.41 respectively, 
consider a range of values for these parameters. Note that this fiducial cross-section 
is, given the recent limits from the CDMS II [21] and XenonlO [22] experiments, an 
order of magnitude smaller than that used in Paper I [TO]. We use the standard values 
for the local circular speed and WIMP density, Vc = 220kms~^ and = 0.3 GeVcm"^ 
respectively. We examine the effect of uncertainties in the local circular speed in Sec. 13.31 
The local WIMP density only affects the amplitude, and not the shape, of the energy 
spectrum. Therefore it only affects the WIMP mass determination indirectly, through 
the number of events detected. We note that the limits from our idealized simulated 
detector are likely to be better than those achievable in reality by a real detector. 

In fig. [2]we plot the 68% and 95% confidence limits on the WIMP mass, m^™, for the 
fiducial detector configuration as a function of the input WIMP mass, m™. Here, and 
throughout, the limits terminate when there is a > 5% probability that an experiment 
will detect no events. Fig. [3] uses the same data, but shows the fractional limits on 

has a lower nuclear recoil acceptance, ~ 30%. 
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Table 1. Dependence of the 95% fractional confidence limits on the WIMP mass, 
(m^™ — m^)/m^, on the energy threshold, Eth, for the benchmark Ge detector, for 
input WIMP masses m^' = 50, 100 and 200 GeV and exposure £ = 3x 10^ kg day. 



Ett (keV) 


50 


(GeV) 
100 


200 





-0.064 , +0.043 


-0.12 , +0.14 


-0.25 , +0.63 


10 


-0.082 , +0.057 


-0.13 , +0.18 


-0.27 , +0.81 


20 


-0.11, +0.073 


-0.16 , +0.22 


-0.31 , +1.1 



the WIMP mass, (m'^™ — m™)/m™. With exposures of £^ = 3 x 10"^ and 3 x 10^ kg day 
it would be possible, with this detector configuration, to measure the mass of a light 
[m^ ~ 0(50 GeV)] WIMP with an accuracy of roughly 25% and 10% respectively. These 
numbers, and the upper limits in particular, increase with increasing WIMP mass, and 
for heavy WIMPs (m^ ^ 100 GeV) even with a large exposure it will only be possible to 
place a lower limit on the mass. For very light WIMPs, < (9(20 GeV), the number 
of events above the detector energy threshold would be too small to allow the mass to 
be measured accurately. 

3.1.1. Energy threshold We now examine the effects of varying the energy threshold, 
Eth (for the fiducial detector configuration and WIMP properties described above). 
Table 1 contains the 95% confidence limits on the fractional deviation of the WIMP 
mass from the input WIMP mass for input WIMP masses of m™ = 50, 100 and 200 GeV, 
energy thresholds Eth = 0, 10 and 20keV and an exposure of £^ = 3 x 10^ kg day. 

As the energy threshold is increased the expected number of events decreases. The 
smaller range of recoil energies also reduces the accuracy with which the characteristic 
scale of the energy spectrum, Eji, and hence the WIMP mass can be determined. The 
effect of varying Eth is smallest for intermediate WIMP masses; for rriy. = 50 GeV (and 
light WIMPs in general) the small Ej^ means that expected number of events decreases 
rapidly as the energy threshold is increased, while for = 200 GeV (and heavy WIMPs 
in general) the large E^, and fiatter energy spectrum, means that the smaller range of 
recoil energies reduces the accuracy with which E-^ can be measured. 

3.1.2. Maximum energy We have previously assumed that recoil events of all energies 
above the threshold energy can be detected. In real experiments there will be a 
maximum energy, -Emax, above which recoils are not detected/analysed. For instance 
for CDMS II [21] -Emax = 100 keV. Fig. H] compares the fractional mass limits for 
-^'max = 100 keV with those found previously assuming no upper limit. The difference 
is very small for light WIMPs [m^ < (9(50 GeV)] increasing with increasing to 
(9(10%) for ~ (9(200 GeV) and £ = 3 x lO^kgday. This refiects the fact that for 
light WIMPs the differential event rate above -Emax is essentially negligible, however this 
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A 

Figure 4. Fractional mass limits as a function of input mass for the fiducial detector 
configuration, with no limit on the energy of recoils which can be detected, (solid lines 
as before) and for a maximum energy i^max = 100 GeV (long dashed). For clarity only 
the 95% confidence limits are displayed in this figure. 

is not the case for heavier WIMPs and finite -Emax reduces the accuracy with which the 
characteristic energy scale of the spectrum can be measured. 

3.1.3. Exposure Fig. [5] shows the fractional hmits on the WIMP mass as a function 
of the exposure, for input WIMP masses of = 50, 100 and 200 GeV, for the 
fiducial detector and dp = 10^^ pb. As the exposure is increased the mass limits (and 
in particular the 95% upper confidence hmit) improve, initially rapidly and then more 
slowly (refiecting the fact that the expected number of events is directly proportional 
to the exposure). An accuracy of ~ ±10% in the determination of the WIMP mass 
can be achieved with an exposure £ = 10^ (10^) kg day for = 50 (100) GeV. For 

= 200 GeV even S = 10^ kg day would not be sufficient to achieve 0(10%) precision. 

3.1.4. Cross-section The expected number of events is directly proportional to both 
the cross-section and exposure. Varying the cross-section is therefore very similar to 
varying the exposure (and hence we do not display a plot of the limits for varying o"p). 
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Figure 5. Fractional mass limits as a function of exposure, f , for the benchmark Ge 
detector and CTp — 10~^pb, for, from top to bottom, m™ = 50, 100 and 200 GeV. 



Unsurprisingly the accuracy with which it will be possible to determine the WIMP 
mass depends sensitively on the underlying cross- sect ion. For instance if dp = 10^^ pb, 
with an exposure of £^ = 3 x 10^ kg day it will be possible to measure the mass with 
an accuracy of ~ ±20%{tS°) = 50 (100) GeV. If dp = lO-^^pb, even for 

= 50 GeV with an exposure of £^ = 3 x 10^ kg day it will only be possible to 
determine the WIMP mass to within a factor of a few and for more massive WIMPs it 
will only be possible to place a lower limit on the mass. 

3.2. Xenon 

We now examine the dependence of the mass limits on the detector target material. In 
fig. Owe compare the fractional limits for the fiducial (Super-CDMS \13\ like) detector 
with Eth = 10 keV with those from a Xe detector with Eth = 4.5 keV (c.f. XenonlO [22] ) 
and lOkeV. The difference in the mass limits for the different detector configurations 
is largest for small exposures, where the number of events observed is small. For light 
WIMPs, as found in Sec. 13.1.11 the energy threshold is important, with the Xe detector 
with Eth = 4.5 keV doing significantly better than both the fiducial Ge detector and also 
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the Xe detector with i?th = 10 keV. For heavier WIMPs, > 100 GeV, in contrast 
to the naive expectation from the dependence of E-^ (see Sec. 12.11) . the Ge detector 
produces shghtly better hmits than the Xe detector with Eth = 4.5 keV. This is because 
of the rapid decrease of the Xe form factor with increasing energy /momentum transfer. 
In fig. [7] we compare, for a Ge detector with Eth = 10 keV and a Xe detector with 
i?th = 4.5 keV, the mass hmits we found before with those which would be obtained if 
the form factor were equal to unity for all energies, F{E) = 1. Without the form factor 
the mass limits are substantially tighter for both Ge and Xe. This is because the event 
rate, in particular at large energies, is increased. The improvement in the accuracy of 
the determination of the characteristic energy, E-^, and hence m^, is greater than if the 
exposure were simply increased so as to increase the expected number of events. For 
instance for Xe, = 200 GeV and S = 3 x 10^ kg day without the form factor the 
one-o" error on the WIMP mass is ~ 3 GeV, while with the form factor included if the 
exposure is increased so as to give the same expected number of events the one-cr error is 
~ 20 GeV. This is because the greater relative abundance of large energy recoils allows 
En to be determined more accurately. With the form factor set to unity the mass 
limits for large WIMP masses are significantly better for Xe than for Ge, as naively 
expected from the WIMP mass dependence of However for any real detector the 
rapid decrease of the Xe form factor with increasing energy /momentum transfer means 
that, contrary to naive expectations, the mass of heavy WIMPs can not be measured 
more accurately with Xe than with Ge (assuming similar threshold energies). This 
conclusion, see also recent discussion by Drees and Shan [H], should also hold for any 
other heavy target. 



3.3. Varying input circular speed, 

Up until now we have assumed that the local circular speed, Vc, is known and equal 
to its standard value of 220kms~^. There is in fact an uncertainty in Vc of order 
±20kms"^ [15] and since E^ depends on both and Vc there is a degeneracy between 
and Vc [10]. Physically, the kinetic energies of the incoming WIMPs depend on their 
mass and velocities. For larger (smaller) Vc the incoming WIMPs have larger (smaller) 
mean kinetic energies than assumed, resulting in larger (smaller) maximum likelihood 
mass values. This can be made more quantitative by differentiating the expression for 
the characteristic energy eq. (IHl): 

— ^ = -[i + KMa)]— ^. (10) 

For an input WIMP mass of m™ = 100 GeV and a 20 km s~^ uncertainty in Vc, this gives 
a ~ 20 GeV shift in the value of the WIMP mass determined. 

In fig. [H]we plot the fractional mass limits for input circular speeds = 200, 220 
and 240 km s"^ for £ = 3 x 10^ kg day. We carry out the likelihood analysis twice, once 
with Vc fixed at 220kms~^ and once with Vc as an additional variable parameter. As 
discussed in Sec. 12.11 in this section we use the fitting function, eq. ([6]), for both the 



12 




Figure 6. Fractional mass limits for the fiducial Ge detector with Eth — lOkeV 
(solid lines for both 68% and 95% confidence limits) and for a Xe detector with 
Eth = 4.5 and lOkcV (dashed and dotted lines respectively) for (from top to bottom) 
£ = 3x 10^, 3 X 10-* and 3 x 10^ kg day. 



input energy spectrum and the likelihood analysis as it is not computationally feasible 
to carry out the full calculation of the energy spectrum for each value of Vc considered 
during the likelihood analysis. 

When the underlying value of is different from the (fixed) value used in 
the likelihood analysis there is, as expected, a significant systematic error in the 
determination of the WIMP mass. For deviations of ±20kms~^ in the underlying 
value of Vc this systematic error increases with increasing m^^ from ~ 10% for small 
to ~ 40% for ^ 200 GeV. The limits are however asymmetric, with the systematic 
error in the upper limits being substantially larger for = 200kms~^. Allowing the 
value of Vc to vary in the likelihood analysis substantially reduces the error, but there 
still appears to be a small systematic shift in the mass limits. 

3.4- Backgrounds 

While future experiments aim to have negligible backgrounds (e.g. Ref. [I3]), non- 
negligible neutron backgrounds would lead to errors in the determination of the WIMP 
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Figure 7. A comparison of the fractional mass limits obtained, for Ge and Xe with 
an exposure £ = 3 x 10^ kg day, with and without the form factor. The solid and long 
dashed lines are for the fiducial Ge detector with Eth = 10 keV and for a Xe detector 
with Eth — 4.5 keV respectively, as in the bottom panel of Fig.[6l The dotted and short 
dashed lines are for the same Ge and Xe detectors, with the form factor artificially set 
to unity, F{E) = 1. 



mass. The size of the errors will depend on the amplitude and shape of the background 
spectrum. In particular if the background spectrum is exponential it can closely mimic 
the shape of a WIMP recoil spectrum (see fig. [9]). 

Motivated by simulations of the neutron background in various dark matter 
detectors [25] we consider two forms for the background 

(i) A flat background energy spectrum from Eth = 10 keV to -Emax = 100 keV , 
parametrised by the total background rate, b,., per tonne year. 

(ii) A exponential energy spectrum from Et^ = 10 keV to -Emax = 100 keV, parametrised 
by the total background rate, b,., per tonne year, and the characteristic background 
energy scale, E\y: 

% In this section we assume that only recoils up to -Emax are detected. 
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Figure 8. Fractional mass limits for the fiducial Ge detector and £ = 3 x 10^ kg day for 
varying circular speed, Vc- In the left panel Vc is fixed as 220 km s^^ in the likelihood 
analysis while in the right panel Vc is allowed to vary in the likelihood analysis (i.e. it is 
an additional fitting parameter). Solid, dotted and dashed fines are for v™ = 220, 200 
and 240kms~^ respectively. 



where 

'dR 



(12) 



The hmits on the WIMP mass for an exposure of S = 3 x 10^ kg day and 
cTp = 10^*^ pb in the presence of a flat background with = 10 and 100 tonne"^ year~^ 
are displayed in fig. [TOl We carry out the likelihood analysis of the WIMP parameters 
twice, firstly neglecting the background and then including the background rate as an 
additional parameter. Neglecting the background leads to a systematic over-estimate of 
the WIMP mass, since the fiat background increases the event rate at large E relative 
to that at small E, so that the best fit energy spectrum has larger E^i, or equivalently 
larger m^. Including the background rate in the likelihood analysis avoids the systematic 
error but, inevitably, leads to larger statistical error in the WIMP mass limits. The 
fractional errors are smallest for ~ (9(50 GeV) and increase for smaller and larger 

+ Note that the backgrounds in real experiments are not expected to be this large. 
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Figure 9. Solid lines are the differential energy spectra including the form factor, 
dR/dE, for (from top to bottom &t E = OkeV) WIMPs with = 50, 100 and 
200 GeV (solid lines). The dotted lines are exponential background spectra with (from 
top to bottom) = 1000 tonne~^year~^ h E\^ — 15keV, &r — 1000 tonne~^year~^ 
& £'b = 25keV and h, = 670 tonne" Vear"! & = 35keV (dotted hues). The 
parameters of the exponential background spectra have been chosen to demonstrate 
that, even when the form factor is included, the WIMP recoil spectra are close to 
exponential and could, in principle, be mimicked by an exponential background. 

WIMP masses. This is because for small WIMP masses the background event rate is 
larger compared with the WIMP event rate, while for large WIMP masses the WIMP 
energy spectrum (exponential with large characteristic energy scale) is closer in shape 
to the flat background spectrum. The systematic (background not included in analysis) 
and additional statistical (background rate included) errors are both at least 10% for 
6r = 10 tonne"^ year~^. For = 100 tonne"^ year~^ the minimum systematic error (for 

~ 50 GeV) is ~ 30%, the shift in the lower limits when the background is included 
in the likelihood analysis is not much larger than for fej. = 10 tonne"^ year~^, however 
the upper limits are increased significantly. 

The limits on the WIMP mass for an exposure of £^ = 3 x 10^ kg day and 
(Tp = 10^^ pb in the presence of a background with an exponential energy spectrum 
with E]j = 20 keV are displayed in fig. [TTl The exponential background spectrum with 
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Figure 10. Fractional mass limits in the presence of a background with a flat energy 
spectrum and total rate br = 10 and 100 tonne" ^ year" ^ (top and bottom panels 
respectively). The dashed (dotted) lines are for when the background event rate is 
not (is) included in the likelihood analysis. The solid lines are the confidence limits 
for zero background. The fiducial Ge detector configuration with E'max = 100 keV, 
£ = 3x10^ kg day and cxp = 10~* pb is used. 



£'b = 20 keV is similar to the WIMP spectrum with ~ 70 GeV. Therefore if the 
background is neglected in the likelihood analysis, for smaller (larger) input WIMP 
masses the WIMP mass is systematically over(under)-estimatcd. For the exponential 
background when the background rate, bj.. and characteristic energy, E'b; arc included in 
the likelihood analysis there arc still large deviations from the zero background limits, 
this is because the WIMP and background spectra can have extremely similar shapes 
and can be degenerate. The fluctuations in the limits with varying reflect the errors 
resulting from this degeneracy rather than a real underlying trend. 

In summary the implications of non-negligible backgrounds for the determination 
of the WIMP mass depend strongly on the shape of the background spectra (as well 
as, obviously, its amplitude). A flat background spectrum will lead to a systematic 
error in the WIMP mass for light WIMPs, which can be avoided, at the cost of larger 
statistical error, by fitting for the background event rate. For heavier WIMPs the fiat 
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Figure 11. As fig. [TO] for an exponential background energy spectrum with 
Eh = 20 keV. 

background is similar in shape to the WIMP spectrum and it is hence more difficult to 
separate the WIMP and background spectra and accurately measure the WIMP mass. 
An exponential background spectrum is similar in shape to the WIMP spectrum and 
would inevitably (even when the background parameters are including in the likelihood 
analysis) lead to increased errors in the determination of the WIMP mass. 

With a single detector it will be difficult to disentangle a WIMP signal (and the 
WIMP mass) from background if the background spectrum has a similar shape to 
the WIMP spectrum (i.e. exponential background, or flat background with a heavy 
WIMP). Multiple targets (for instance Ge and Si as used by CDMS II [21]) would help 
due to the dependence of the WIMP spectrum on the mass of the target nuclei. See 
ref. [26] for Monte Carlo simulations using CaW04 and ZnW04. Detectors composed of 
very different targets (e.g. Ge and Xe) would likely have different background spectra 
however. 
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4. Summary 

We have studied how the accuracy with which the WIMP mass could be determined by a 
single direct detection experiment depends on the detector configuration and the WIMP 
properties. Specifically, we investigated the effects of varying the underlying WIMP mass 
and cross-section, the detector target nucleus, exposure, energy threshold and maximum 
energy, the local circular speed and the background event rate and spectrum. 

The accuracy of the mass limits is most strongly dependent on the underlying 
WIMP mass and the number of events detected. For light WIMPs (mass significantly 
less than that of the target nuclei) small variations in the WIMP mass lead to significant 
changes in the energy spectrum. Conversely for heavy WIMPs the energy spectrum 
depends only weakly on the WIMP mass. Consequently it will be far easier to measure 
the WIMP mass if it is light than if it is heavy. The number of events detected is directly 
proportional to both the exposure and the cross-section, therefore these quantities have 
the greatest bearing on the accuracy of the WIMP mass determination. For our baseline, 
SuperCDMS [13] like, Ge detector with negligible background and energy threshold 
Eth = lOkeV for a WIMP-proton cross-section of cxp = 10~^pb, a factor of a few below 
the current exclusion limits from the CDMS II [21] and XenonlO [22] collaborations, 
with exposures of £^ = 3 x 10"^ and 3 x 10^ kg day it would be possible to measure 
the mass of a light [m^ ~ (9(50 GeV)] WIMP with an accuracy of roughly 25% and 
10% respectively. These numbers, and the upper limits in particular, increase with 
increasing WIMP mass, and for heavy WIMPs {m^ ^ 100 GeV) even with a large 
exposure it will only be possible to place a lower limit on the mass. For very light 
WIMPs, < (9(20 GeV), the number of events above the detector energy threshold 
would be too small to allow the mass to be measured accurately. If dp = 10~^ pb, with an 
exposure of £ = 3 x 10^ kg day it will be possible to measure the mass with an accuracy 
of ~ ±20% itS) for m}^ = 50 (100) GeV. If cXp = 10"^° pb, for = 50 GeV, even 
with an exposure of £ = 3 x 10^ kg day it will only be possible to determine the WIMP 
mass to within a factor of a few and for more massive WIMPs it will only be possible 
to place a lower limit on the mass. 

The energy threshold, E'th, and the maximum energy, -Emax, above which recoils 
are not detected/analysed also affect the accuracy with which the WIMP mass can be 
determined. Increasing E^-^ (or decreasing -Emax) not only reduces the number of events 
detected, but also reduces the range of recoil energies and the accuracy with which the 
characteristic energy of the energy spectrum, Eji, and hence the WIMP mass, can be 
measured. The effect of increasing Eth is smallest for intermediate WIMP masses. For 
light WIMPs the small E^ means that the expected number of events decreases rapidly 
as the energy threshold is increased, while for heavy WIMPs the large Eji, and flatter 
energy spectrum, means that the smaller range of recoil energies reduces the accuracy 
with which Eji can be measured. The effect of reducing the maximum energy (from 
infinity to 100 keV) is very small for light WIMPs as the differential event rate above 
Ema.x = 100 keV is negligible, however for heavy WIMPs the fractional mass limits can 
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change by 0(10%). 

The relative capabihties of different detectors to determine the WIMP mass depend 
not only on the WIMP and target masses, but also on their energy thresholds. The 
WIMP and target mass dependence of the characteristic energy scale of the recoil 
spectrum suggests that heavy targets will be able to measure the mass of a heavy 
WIMP more accurately, however the rapid decrease of the nuclear form factor with 
increasing momentum transfer which occurs for heavy nuclei means that this is in fact 
not the case. 

If the WIMP distribution on the ultra-local scales probed by direct detection 
experiments is smooth, then the uncertainties in the detailed shape of the local velocity 
distribution lead to relatively small changes in the shape of the differential event 
rate [I7l|18], and hence a relatively small, [C(5%)], systematic uncertainty in the WIMP 
mass [1^. There is however an uncertainty in the local circular speed, fc, (and hence 
the typical speed of the WIMPs) of order ±20kms~^ [15] and since E^i depends on 
both and Vc this leads to a degeneracy between and Vc [ID]. For deviations of 
±20kms^^ in the underlying value of Vc this systematic error increases with increasing 
m'^ from ~ 10% for small m'^ to ~ 40% for ^ 200 GeV. 

The assumption of a smooth WIMP distribution may not be valid on the sub milli- 
pc scales probed by direct detection experiments (see discussion in Paper I [10]). If 
the ultra-local WIMP distribution consists of a finite number of streams (with a priori 
unknown velocities) then the recoil spectrum will consist of a number of (sloping due 
to the energy dependence of the form factor) steps. The positions of the steps will 
depend on the stream velocities, the target mass and the WIMP mass. In this case 
multiple targets would be needed to extract any information on the WIMP mass. Drees 
and Shan [11] have recently demonstrated that with multiple targets it is in principle 
possible to constrain the WIMP mass without making any assumptions about the WIMP 
velocity distribution. 

Future experiments aim to have negligible backgrounds, however, non-negligible 
neutron backgrounds would lead to errors in the determination of the WIMP mass. 
The size of the errors will depend on the amplitude and shape of the background 
spectrum. If the background rate is not negligible compared with the WIMP event 
rate it will be difficult to disentangle a WIMP signal (and the WIMP mass) from the 
background if the background spectrum has a similar shape to the WIMP spectrum (i.e. 
exponential background, or flat background with a heavy WIMP). The uncertainties 
from backgrounds could be mitigated by using multiple targets (see e.g. Ref. [26]). 
however detectors composed of very different targets (such as Ge and Xe) would be 
unlikely to have the same background spectra. 
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